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The above-noted amendment of claims 1-21 as well as 
the submission of an Abstract and revisions to the 
Specification, are respectfully submitted prior to initiation of 
the prosecution of this application in the U.S. Patent and 
Trademark Office. 

The above-noted new claims are respectfully submitted 
in order to more clearly and appropriately claim the subject 
matter which applicant considers to constitute his inventive 
contribution. No new matter is included in these amendments. 
In addition, the revisions to the Abstract and Specification are 
submitted in order to clarify and correct the Abstract and 
Specification and to conform them to all of the requirements of 
U.S. practice. No new matter is included in these amendments. 

In view of the above, it is respectfully requested 
that these amendments now be entered, and that prosecution on 
the merits of this application now be initiated. If, however, 
for any reason the Examiner does not believe such action can be 
taken, it is respectfully requested that the Examiner telephone 
applicant's attorney at (908) 654-5000 in order to overcome any 
objections which the Examiner may have. 
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METHOD AND DEVICE APPARATUS FOR ASSIGNING 

WEIGHTING COEFFICIENTS FOR PERFORMING 
ATTITUDE CALCULATIONS WITH A STAR SENSOR 

General technical domain 

CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] The present application is a national phase entry 

under 35 U.S.C. § 371 of International Application 

No. PCT/FR04/00675 filed March 19,- 2004, published in France, 

which claims priority from French Application No. 03-03407, 

Publication No. 2852687 filed March 20, 2003, all of which are 

incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

[0002] This invention relates to the domain of methods for 

determining the orientation of objects in space. 

[0003] More precisely, it relates to a method of assigning 

weighting coefficients to measurements of a succession of 

stars acquired by a star sensor in order to determine a 

spatial orientation. 

State of the art 

[0004] An increasing number of artificial satellites are 
orbiting around the Earth or are launched into space. 
Obviously, it is important to know the position of satellites 
in space, but it is also essential to know their attitude, in 
other words their orientation with respect to an inertial 
coordinate system in the celestial dome. 

[0005] A star sensor is an instrument that supplies its own 
attitude measurement calculated from direction measurements of 
stars that it carries out, or that directly transmits the list 
of measurements of elementary directions to a client device, 
for example such as the central software on a satellite on 
which the sensor is located. 

[0006] Figure 1 shows a satellite 1, for example in an 

orbit 2 around a planet 3, for example the Earth, at an 
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altitude 13 above the planet surface. The satellite 1 
comprises a star sensor 4 on one of its walls facing the stars 
in the celestial dome, and with a view field diagrammatically 
shown by the trunk of the cone 5. The circle 8 represents the 
part of the celestial dome observable in the view field 5. 
The stars 10 and 11 diagrammatically represent some of the 
stars that can be observed by the sensor 4 . It can be 
understood that a plurality or a large number of stars can be 
observed at the same time within the circle 8. The stars 9 
and 12 represent stars in the celestial dome outside the view 
field of the sensor. 

[0007] The satellite 1 displaces in the orbit 2 along the 
arrow 6. Consequently, it can be understood that the circle 8 
also displaces on the celestial dome along the arrow 7. Thus, 
the star 9 used to be in the circle 8, and the star 12 will be 
in the circle 8 in the future of the satellite displacement. 
[0008] Most frequently, the attitude of the sensor 4 or the 
satellite 1 is calculated using only some of the stars present 
in the view field 8 of the star sensor 4. For example, about 
ten stars may be chosen to make an attitude calculation. The 
direction of the star 10 or 11 with respect to a coordinate 
system 15 related to the sensor 4 is measured in the attitude 
calculation. Since the direction of the star 10 or 11 with 
respect to a known inertial coordinate system 14 is also 
known, the direction of the coordinate system 15 (and 
therefore of the satellite 1 or the sensor 4) with respect to 
the coordinate system 14 of the celestial dome may be deduced. 
[0009] The precision of the attitude estimate depends on 
the choice of stars to be measured for the direction 
calculation . 

[0010] The choice of stars is also important when star 
direction measurements are transmitted directly and are used 
by the client device of the star sensor. 
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[0011] In some cases, the stars may be associated with a 

weight instead of being selected or rejected, which weights 
the importance to be assigned to them in subsequent 
processing . 

[0012] In other cases, the star sensor itself makes a first 

selection of stars, makes measurements on these selected stars 
and may send them accompanied by a weight or transformed into 
a global attitude estimate. 

[0013] The measurements may also originate from different 
star sensors onboard the same satellite. 

[0014] In all cases, a selection is a set of potential or 

previously made measurements of star directions, for which the 
validity date is identical. As we have already said, weights 
may be associated with each measurement. Over time, a star 
sensor 4 generates a succession of selections, measurements 
being made on all or some of the selected stars. 
[0015] Traditionally, stars chosen as being the best stars 

will be selected for as long as possible. In the early dates 
of stars sighting, it was difficult enough simply to lock onto 
a star. Therefore, once it was found it was not released. 
[0016] One of the criteria for choosing stars in a 
selection is particularly the magnitude of each star that can 
be observed in the view field. Thus, the brightest stars 
should be chosen, namely stars with the lowest magnitude. One 
of the other selection criteria is the distance of the star 
from the optical axis of the view field. More precisely, 
Figure 2 shows the circle 8 moving along the arrow 7 across 
the celestial dome. Choosing stars remote from the optical 
centre 83 and located within a reference area 81 within the 
circle 8 can improve the precision of the estimate depending 
on the optical axis 83 of the sensor. The axis 83 is the most 
likely to be affected by measurement errors. Stars located in 
an area reference 82 will be close to the optical axis 83 at 
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one time or another and will provoke measurement errors. A 
star that is too close to the optical centre 83 will be 
re j ec ted and another star further from the optical centre 83 
will then be selected. 

[0017] At the moment, selection methods give excessive 
importance to these two criteria, with the result that stars 
are not changed frequently. 

[0018] However, the above methods have disadvantages. 
[0019] Some stars, for which the direction is particularly 

badly estimated by sensor 4, disturb the global attitude 
estimate throughout the duration in which they are selected. 
This duration may be fairly long when the attitude of the star 
sensor 4 varies only slightly and when its view field is 
fairly wide. Figure 3 diagrammatically shows such a 
situation. Figure 3 shows the error on the attitude of the 
sensor 4 as a function of time. Changes 30 in the error level 
of the graph are due to star selection changes for 
determination of the attitude. It is seen that errors due to 
the choice of a selection are relatively long compared with 
the oscillations 31 on each plateau of the curve, which are 
due to observation errors of each star in the selection. It 
is observed that the errors 31 cancel out due to averaging 
over a relatively short time compared with the time during 
which each selection is observed. In other words, the noise 
due to oscillations 31 can easily be filtered by existing 
techniques for processing of data received from star sensors 
by client devices, since it is within the high frequencies of 
the frequency spectrum of the signal acquired by the sensor. 
The fact that the selection changes only infrequently results 
in low frequency noise that is difficult to filter. 
[0020] More generally, methods of assigning a weighting 
coefficient according to prior art are incapable of 
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controlling the noise dispersion with time due to each 
selection . 

[0021] Furthermore, some phenomena (for example distortion) 
vary depending on the position of stars in each selection in 
the view field. 

[0022] More generally, methods for assigning weighting 
coefficients according to prior art are incapable of 
controlling these space-time dispersion phenomena for stars to 
which high weights are assigned. 

Presentation of the invcntion SUMMARY OF THE INVENTION 
[0023] The invention proposes to overcome these 
disadvantages . 

[0024] One of the purposes of the invention is to divulge a 
method for assignment of weighting or weighting coefficients 
for a calculation of a spatial orientation in order to control 
measurement errors of a star sensor. 

[0025] Another purpose of the invention is to vary the 
dispersion of stars for which the direction measurements are 
used either by the star sensor itself or by the client device 
of the star sensor. 

[0026] Another purpose of the invention is to propose a 
method of enabling space-time dispersion of star selections. 
[0027] Finally, another purpose of the invention is to 

propose a method that is capable of taking account of 
previously made selections in order to spread low frequency 
noise related to star selections more or less throughout the 
spectrum. 

[0028] To achieve this, the invention proposes a method for 

assignment of weighting coefficients to measurements of a 
succession of stars acquired by a star sensor related to a 
client device to determine a spatial orientation, 
characterised in that higher or lower preference is given to 
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refreshment of the positions of measurements with the highest 
weights and / or stars on which these measurements are made by 
the star sensor and / or its client device, so as to displace 
part of the power of the error associated with the set of star 
measurements within the frequency spectrum. 
[0029] This invention also advantageously includes the 
following characteristics, taken alone or in any technically 
possible combination: in the calculation of the weights of 
measurements in a current selection, the reinforcement or 
attenuation takes place as a result of applying a distance 
weight associated with each measurement in the current 
selection and characteristic of an average distance between 
firstly the said measurement and secondly the measurements in 
the previous selections and the other measurements in the 
current selection; the distance weight associated with the 
current selection measurement is calculated as a weighted 
average of the corresponding distances between firstly the 
said measurement, and secondly the previous selection 
measurements and the other measurements in the current 
selection respectively; the weighting coefficient associated 
with the distance between a first measurement in the current 
selection and a second measurement in a previous selection or 
another measurement in the current selection includes a memory 
coefficient associated with the said second measurement, and / 
e^ -and/or the weight of the second measurement if it belongs 
to a previous selection or a temporary weight if it belongs to 
the current selection; the distance calculation combines the 
angular distance between the two measurements, and an identity 
distance that depends on the difference in nature of the two 
stars for which the measurements are being made; the memory 
coefficient of a measurement m± at a time t is defined using 
the following formula: 

Mem (mi/ 1) = \ixU' [t ~ T<m i )] f where 
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• T(m k ) is a validity date of a measurement m k 

• p and n are constants. 

[0030] a-A charge is assigned to each star for which a 
measurement is made, the charge summarising the weights 
assigned to the measurements made on the said star in the 
past, attenuated by the passage of time; 

the charge of the star e p is defined at an instant T by the 
following formula : 

Cha (e p , T) = jc MemCnij/T)] 

/=/>+! 

where Mem(rrii/T) is the memory coefficient of the measurement mi 
at time T, E(mi) is the star on which the measurement m± is 
made, and A (m±) is the measurement weight m±; 

the charge associated with a star to which a measurement in 
the current selection is related is updated before it is used 
in the calculation of the weight associated with a 
measurement, using a coefficient that depends on the 
difference A between the current date and the last update date 
for this charge; 

the coefficient may be a factor and is in the form n~ A , where 
n is a constant; 

the coefficient may be additive and is in the form - pxA , where p 
is a constant; 

after calculating the weight associated with a measurement in 
the current selection, the charge associated with the star for 
which this measurement was made is updated; 

the update is made by adding the weight associated with the 
measurement ; 

a random function is used in the calculation of the weights; 
the calculation of the distance weight is iterated with a 
temporary weight for measurements in the current selection, 
the distance weight being used to calculate a new weight 
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itself used to calculate a new distance weight and so on, 
until convergence towards a final weight; 

the digital values of the method are saved in memory and 
processing means of the sensor and / or the client device; 
the renewal rate of stars with a large weight is increased by 
increasing the frequency of measurements of the star sensor 
and / or the client device; 

the dispersion of the complete new selection is used directly 
in the weights, using processing means related to the sensor 
and / or client device; 

processing means related to the sensor and / or the client 
device comprising a neurone structure are used to directly 
affect dispersion in the weights. 

[0031] The invention also relates to the device for 
implementing this method. 

Presentation of the figures 

[0032] Other characteristics, purposes and advantages of 
the invention will become clear from the following description 
that is purely illustrative and is in no way limitative, and 
that should be read with reference to the attached drawings in 
which : 

Figure 1, already commented upon, diagrammatically shows an 
artificial satellite in orbit about a planet, and comprising a 
star sensor; 

Figure 2, already commented upon, diagrammatically shows a 
view field of a star sensor on the celestial dome; 
Figure 3, already commented upon, diagrammatically shows a 
graph of the error on the attitude as a function of time in 
methods for assignment of weighting coefficients according to 
the state of the art; 
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Figure 4 diagrammatically shows a graph of the error on the 
attitude as a function of time in a method for the assignment 
of weighting coefficients according to the invention. 

Detailed dcocription DETAILED DESCRIPTION OF THE INVENTION 
[0033] In some of its calculation phases, a star sensor 4 

shown in Figure 1 chooses to select specific stars (10 and 11 
on Figure 1) that it detects within its view field 8, either 
to use it itself in the calculation of the attitude, or to 
transmit the measurement of its direction directly to the 
client device, for example the satellite 1. 
[0034] The selection then involves the assignment of a 
weight to each star 10 and 11 that characterises the potential 
of the star with respect to the needs of the client device 1 
of the star sensor 4. 

[0035] According to other methods, the stars are used in 
the calculation of the attitude with weighting, or a weight is 
transmitted to the client device 1 with the measurement of the 
star direction. All stars may then be transmitted, and a star 
transmitted with a zero weight is equivalent to an unselected 
star . 

[0036] In any case, a weight is associated with a 

measurement of a star direction at one time or another. This 
weight characterises the potential of the measurement to 
satisfy the needs of the client device 1 of the star sensor 4. 

[0037] We will now describe a first possible method for 
assigning weighting coefficients to measurements of a 
succession of stars acquired with a star sensor in order to 
determine a spatial orientation. According to this first 
method, the weights are assigned to the measurements 
individually . 
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[0038] According to a first method, the importance assigned 
to measurements is controlled as a function of whether or not 
they relate to new stars, which may be called the time 
dispersion of star measurements. It takes account of previous 
selections. The dispersion effect of measurements is 
amplified . 

[0039] According to a first possible embodiment, the 
consecutive duration during which a particular star is 
selected is reduced during observation of the stars. The 
effects of such a situation are shown diagrammatically in 
Figure 4. Thus, the error on the attitude of the sensor 4 in 
Figure 1 as a function of time is plotted in Figure 4. As in 
Figure 3, the changes 30 in the error level on the graph are 
due to changes in star selections made to determine the 
attitude. It is found that the errors due to the choice of a 
selection are relatively shorter than the errors in Figure 3. 
They are of the same order of magnitude as the errors 31 due 
to observation errors of each star in the selection. In terms 
of the spectrum, a higher frequency of changing the star 
selection is equivalent to replacing low frequency noise 
(induced by stars made noisy by oscillations 31 and kept for a 
long time) by higher frequency noise (induced by a more 
frequent change of stars used in the attitude estimate) . More 
precisely, this is equivalent to spreading the power at very 
low frequency within a wider frequency range. This operation 
may be particularly interesting if the client device of the 
star sensor has higher frequency measurement capabilities 
(which is the case when it uses high performance gyroscopes) . 
Noise can be easily filtered. 

[0040] Conversely, it can be understood that if high 
frequency filtering capabilities of the client device are 
weak, and if the high frequency noise is bothersome, it may be 
useful to give preference to low frequencies, by keeping each 
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star for a longer period. Therefore, there will be a mode for 
using the method controlling a different dispersion of 
measurements in time. 

[0041] According to another possible method according to 
the invention, the dispersion of star measurements in space 
will be controlled. 

[0042] According to a first possible embodiment, when stars 
are being observed, the selected stars are geometrically 
dispersed in the view field. This controls measurement 
distortion. This dispersion relates to a selection made at a 
given instant . 

[0043] It can also be understood that according to another 

embodiment, it is required to not disperse the measurements in 
space . 

[0044] Space - time dispersion of measurements will be 
combined to control the error due to each selection in the 
orientation calculation . 

[0045] Consequently, in calculating the weight associated 
with a measurement, the invention uses the average distance of 
the corresponding measurement with measurements in previous 
selections and the present selection. If it is required that 
the measurements to which a high weight is assigned should 
last for a long time, then an attempt will be made to minimise 
this average distance. On the other hand, an attempt will be 
made to maximise the distance if it is required that the 
measurements should be renewed frequently. 

Definitions 

[0046] A "distance weight" means an average distance of a 
measurement with measurements in previous and present 
selections. The distance weight is then used in the 
calculation of the weight, which takes account of other 
considerations (for example the position of the measurement on 
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the view field or the magnitude of the star considered) . One 
possible expression of the distance weight is given in the 
remainder of this description. The first step is to define 
some of its possible variables. 

[0047] Let {mi,m 2 , ... m N } be the series of measurements of 
previous and present star directions. Two distinct elements 
in this sequence may correspond to the same physical star 
measured at two different instants. 

[0048] The sub-series {mi,m 2 , ... m p }, where P < N, is the 

series of stars in the present selection. 

[0049] S(mi,mj) characterises the difference between 

measurements m± and m 3 and uses two terms with very different 

natures, for example by their product: 

DA (mi, m-j) - this is the angular difference in the 
coordinate system of the star sensor 4 (or any other 
corresponding distance, for example the geometric 
distance on the detection matrix of the same sensor) . 
Characterising this angular difference provides a means 
of dispersing or grouping measurements that will be 
assigned a high weight on the sensor matrix. For 
example, this provides a means of compensating for some 
optical distortion errors. 

DE (mi, m-j) - this is the identity difference, in other 
words the difference in characteristics of stars for 
which two measurements are made. For example, this 
difference may be equal to 0 if the two stars are 
identical, and otherwise to a positive value (two 
measurements corresponding to identical stars can give 
different values of the distance DA, due to movement of 
the star sensor along the direction 6 in Figure 1) . 
Characterising this distance makes it possible to replace 
stars for which measurements are made more or less 
frequently. This also provides a means of giving 
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preference to stars that have spectra or magnitudes that 
are different from (or similar to) previously selected 
stars . 

[0050] One possibility consists of using the identity 
difference to characterise the difference between the current 
measurement and measurements made in previous selections, and 
using the angular difference to characterise differences 
between measurements in the current selection. 
[0051] The distance may possibly be a function of one or 
both of the previous distances. If the function is decreasing 
both as a function of DA and DE, a proximity weight will be 
calculated instead of a distance weight. This will then be 
used in the formula for calculating the weight instead of. or 
in addition to the distance weight. 

[0052] For example, we will choose 8(m±,mj) = qx (e DA(m j ' m i ) - 

1) xDE(m-j,mi) , where q is a constant, DA(m if mj) is the angular 
distance between the two measurements, and DE(mi,mj) is equal 
to 0 if the stars for which the two measurements were made are 
identical, and otherwise 1. Therefore, S(m ± ,nij) is not zero 
only if the stars are different. 

[0053] The previously calculated weight A (mi) , which in 
particular uses the distance weight, is associated with each 
measurement m± in previous selections. 

[0054] The weight of previous measurements used in the 
calculation of the distance weight corresponding to a 
measurement irij of the current selection is reduced, because 
time has erased the memory of this measurement. The result of 
this reduction is the remanent weight. 

[0055] The remanent weight may be calculated as being the 
product of the weight and a memory coefficient Mem(m±/t) that 
characterises the stored memory of the measurement itu at date 
t. The measurement memory coefficient iru for the measurement 
rrij at date T(mj) is written as Mem (m ± /T (rrij) ) . 
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[0056] The memory coefficient Mem (m±/t) decreases when the 
measurement mi becomes more remote in time. This reflects the 
gradual erasure of star measurements from memory. 
[0057] For example, the following product will be used: 
Mem(m±/t) = iixir It ~ rfl V 1 ; where 

• T(m k ) is the validity date of measurement m k 

• p and n are constants. 

[0058] The values of constants may depend on the difference 
S(mi,mj) . Thus, by varying p, the memory coefficient may 
change to the value 0 when [t-T (m±)] exceeds a certain value. 
This erases the memory of the measurement. The coefficient 
may also change to the value 0 when [t-T (m±)] is less than a 
certain value. This provides a means of forcing the selection 
of stars present in very recent selections. It also provides 
a means of excluding the distance weight of measurements in 
the current selection from the formula. 



Calculation of the distance weight 

[0059] The distance weight associated with the measurement 
irij is an average over the set of measurements in previous 
selections and measurements in the current selection (except 
m-j) . 

[0060] The average may be weighted by remanent weights of 

measurements used in the average. 

[0061] For example, the average will be calculated in the 
following form: 



P(mj) = 



2 [Aim, )xMem{m i I T{m j ))f xSfaj , m, ) * 



£ [Aim, )xMem{m l I T(mj ))f 



1=1 



\lq> 



where 



and 9 are constants. If Ci and q> are equal to 1, then it 
is a linear average of the distances. If q> is equal to 



14 



Application No .: 10/549, 737 



REGIM 3.3-065 



2, it is a root mean square of the same distances (in 
principle Q is equal to 1) . If cp < 1, the separation 
between pairs of very close measurement is accentuated. 
[0062] If the measurement irij is included in the present 
selection, the weight A(mj) will not use the distance weight. 
[0063] In practice, the weight of a measurement in the 

current selection may be calculated from an arbitrary distance 
weight. An iterative method may be used after this first 
calculation: weights associated with these measurements in 
the current selection are used in the calculation of distance 
weights for these measurements, and take account of the 
distance weight calculated in the previous iteration. All 
distance weights of measurements in the current selection are 
thus recalculated until they converge. 

[0064] Another possibility consists of not using the other 
measurements in the current selection in the calculation of 
the distance weights of the current selection. This is 
equivalent to choosing Mem(m±/t) = 0 for all values of i such 
that mi is a measurement in the current selection. 
[0065] The calculation of the distance weight can give a 
weight for each star measurement. 

[0066] For example, the weight may be calculated as 
A(nij) =e~ M * Vfnij) * P (nij) , where M is the magnitude of the star for 
which the measurement irij is made, ^(mj) is the angular 
distance between the measurement irij and the optical axis of 
the star sensor, and P(mj) is the distance weight associated 
with the measurement irij . 

[0067] Thus, the stars with the greatest weights are given 
preference in a selection. 

[0068] It can then be understood that the space - time 
dispersion of selections is controlled by varying the 
different distance weight variables, which results in a 
different final weight to the measurements. The selections 
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may be change more or less quickly, and the selections are 
more or less dispersed among themselves on the matrix of the 
sensor, as a function of the parameters of the chosen 
variables . 

[0069] We will now describe a variant of the method in 
which the distance weight of the measurements in previous 
selections is calculated recursively, to limit the amount of 
information to be memorised and to simplify the calculations. 
[0070] In this variant of the method, the identity distance 

is used to define proximity between measurements. 

Dcf initiono 

[0071] Depending on the variant of the method, the 

information that summarises the past characterises the stars 
on which the measurements in previous selections were made, 
rather than the measurements themselves directly. This 
information can be summarised as a pair C p associated with 
each star e p , for which the measurement formed part of a 
previous selection defined by: 

C p = {Cha, date}, where 
date is the validity date of the charge Cha, which varies with 
time. The charge Cha is the accumulation of remanent weights 
associated with the set of measurements at the instant 
considered, that in previous selections applied to star e p . 
Therefore the charge associated with star e p at time T is 
theoretically defined as follows, using the same notations as 
above and if E(m±) is the star to which the measurement itu 
applied: 

N 

Cha (e p , T) = j [Aim, ) x MemOn/T)] 

E(m f )=e p 

[0072] The term inside the sum can be multiplied by a 
function of the angular distance between the current 
measurement applied to star e p and the measurement mi, to take 
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account of the spatial dispersion between successive 
selections (the target area in the sky possibly changing under 
the effect of the movement of the satellite) . For example, 
this function may be in the following form: 

-r}xd(m(e ) ,m ) 
" pi 

where m(e p ) is the measurement in the current selection 
applicable to star p, d is the distance between the two 
measurements in the argument and r| is a coefficient. 
[0073] Instead of a charge, a memory associated with the 
star may be calculated equal to the average of the memory 
coefficients associated with the measurements that applied to 
the star considered weighted by their weights, for example: 



^[A(m f )xMem(m g /T)] 



i=\ 

SOU(e p ,T) = ' 



1=1 

£(m f .)=e p 

[0074] This memory is a "normalised charge" that will be 
used in the calculation of the weight. 

Propagation of charges 

[0075] Charges are propagated as follows at the time that 
the measurements in the current selection are processed. 
[0076] Charges are refreshed and associated dates are 
updated. Charges associated with stars that form part of 
previous selections are diminished to translate erasure from 
memory. This update may possibly concern only stars in the 
current selection, but all information related to other stars 
must be kept. Charges associated with stars that do not form 
part of any previous selection remain zero. Dates associated 
with charges of stars in previous selections change to the 
date of the current selection. 
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[0077] This reduction in the charge may be made if A is the 
time difference defined as follows: 

A = current selection date - date previously associated 

with the charge, 
by using one of the following two methods: 

1. multiplication of a previous charge by a factor: 

n" A , where n is a constant 

2. addition of a factor to the previous charge: 

-pxAr where p is a constant. 
[0078] The charge thus obtained is forced to 0 if the 
result of the previous operation is negative. 

[0079] The charge is calculated without approximation when 
the previous method has been chosen and the memory coefficient 
is also in exponential form. In other cases, the calculation 
is approximate. 

[0080] The next step is to assign distance weights. The 
distance weight assigned to each measurement in the current 
selection is equal to the refreshed charge associated with the 
star for which the measurement considered is made. If no 
previous measurements were made for the star, its initial 
charge is zero and therefore its distance weight will also be 
zero . 

[0081] The next step is to calculate the weights. The 
distance weights calculated in the previous step are used to 
calculate weights associated with measurements in the current 
selection, for example using the formula A (uij) =e~ M * if/ (itij) x P (itij) . 
[0082] The charges are then updated. Charges of stars to 
which measurements in the current selection are made are 
increased by the value of the corresponding measurement 
weights. For each star e h that was not previously 
characterised by a pair C h = {Cha, date}, the date is fixed at 
the current selection date. For all stars concerned, the date 
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may be assigned in this step rather than in the charge 
refreshment step. 

[0083] Finally, charges that are too small are eliminated. 

This step is optional. Charges below a minimum value may be 
eliminated, which is equivalent to forgetting that at least 
one selection was made on the star in the past. 
[0100] A random function, for example a random Gaussian or 
uniform variable may be used for calculating the weights. 
This contributes to space - time dispersion of stars, in this 
case by chance. 

Uoing the method 

[0101] The distance weight or the charge pair is calculated 
assuming that information is memorised about measurements in 
previous selections, for example chosen from the following for 
each measurement: 

• The measurement date, 

• The weight finally assigned to the measurement, 

• The direction of the measurement, identified by its 
angle in a coordinate system 15 related to the sensor 4, 
or by vector in a coordinate system related to its 
matrix, 

• Some characteristics of the star used to make the 
measurement . 

[0102] If the recursive method is used, the pair C p = {Cha, 

date} is kept for each star for which a set of measurements 
was made for which the memory was not erased. 
[0103] Therefore, memory and processing means are 
associated with the sensor and / or and/or its client device, 
namely the satellite. 

[0104] If the available memory limit is reached, 

information about measurements for which the memory is weakest 
can be deleted. The criterion can apply to: 
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• the product of the memory associated with the 
measurement by the weight of the same measurement; 

• the charge, particularly if the recursive method is 
used. 

[0105] We will now describe a second possible method for 
assigning weighting coefficients to measurements of a 
succession of stars acquired by a star sensor in order to 
determine a spatial orientation. According to this second 
method, the best selection of stars for determining the 
spatial orientation will be chosen directly. We will no 
longer consider stars individually to determine which have the 
greatest coefficients before grouping them in selections. The 
dispersion of a complete new selection will be used directly 
in the weights, using processing means connected to the 
renewal rate of stars with a large weight is increased by 
increasing the frequency of measurements of the star sensor 
and / or the client device; 

[0106] According to this method, the processing means 
connected to the sensor may comprise a neurone — otructurc neural 
network . The neurone stru ct ure neural network, can assign 
weights to stars to determine the orientation of the sensor or 
the satellite, for example after a learning process. The 
learning criterion on the neurone network neural network can 
use distance coefficients, average distances, charges and— /—or 
memory defined above. It may also use the average of distance 
weights, possibly weighted by the weights of stars (or by 
temporary weights that do not take account of distance 
weights) for the stars in a selection. 

[0107] In both methods, the selected star dispersion 

phenomenon may be amplified by forcing the star sensor to 
operate at a frequency higher than is strictly necessary. 
Provided that the previous methods are used, the number of 
selections used will be increased, causing additional 
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spreading of noise on the spectrum. The advantages thus 
obtained could compensate the increase in noise induced by 
shortening of the integration of light information by the star 
sensor, from the point of view of the client device in the 
star sensor. 

[0108] Star measurements can then be filtered and then sub- 
sampled before use, to return to the required frequency. In 
principle, filtering involves the calculation of the attitude 
quaternion. The method is efficient even if the previous and 
proximity weights are not used in the criterion. 
[0109] The above developments are equally applicable 
firstly to a star set detected and selected regularly (for 
example once every second) , at the same time (as is the case 
for conventional operation of a star sensor with a CCD 
detector) , or secondly to a star set in which the stars are 
detected irregularly (use of new types of detectors) and in 
which selections are produced irregularly (as is the case for 
APS (Active Pixel Sensor) star sensors. 

[0110] The methods may be used at several levels, possibly 

simultaneously, for example the selection level of the stars 
on which the stellar sensor (s) will make their measurements; 
weighting of measurements made at the star sensor (s); 
weighting of measurements made at the measurements client 
device . 
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